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ABSTRACT 

We use a high-resolution TV-body simulation to study how the formation of 
cold dark matter (CDM) halos is affected by their environments, and how such 
environmental effects produce the age-dependence of halo clustering observed 
in recent TV-body simulations. We estimate, for each halo selected at redshift 
z = 0, an 'initial' mass M defined to be the mass enclosed by the largest sphere 
which contains the initial barycenter of the halo particles and within which 
the mean linear density is equal to the critical value for spherical collapse 
at z = 0. For halos of a given final mass, Mh, the ratio Mj/Mh has large 
scatter, and the scatter is larger for halos of lower final masses. Halos that 
form earlier on average have larger Mj/Mh, and so correspond to higher peaks 
in the initial density field than their final masses imply. Old halos are more 
strongly clustered than younger ones of the same mass because their initial 
masses are larger. The age- dependence of clustering for low-mass halos is 
entirely due to the difference in the initial/final mass ratio. Low-mass old 
halos are almost always located in the vicinity of big structures, and their 
old ages are largely due to the fact that their mass accretions are suppressed 
by the hot environments produced by the tidal fields of the larger structure. 
The age-dependence of clustering is weaker for more massive halos because 
the heating by large-scale tidal fields is less important. 
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1 INTRODUCTION 

In the standard cold dark matter (CDM) paradigm of structure formation, virialized CDM 
halos are both the building blocks of the large-scale structure of the Universe and the hosts 
within which galaxies are supposed to form. During the last decade, the properties of CDM 
halos, such as their internal structure, formation histories and clustering properties, have 
been studied in great detail using both iV-body numerical simulations and analytical models. 
These studies have shown that the halo bias depends strongly on halo mass, in the sense 
that more massive halos are more strongly clustered (e.g., Mo & White 1996; Jing 1998; 
Seljak & Warren 2004). This mass dependence of the halo bias has played a crucial role 
in understanding the correlation function of both dark matter and galaxies, via the halo 
model (e.g., Mo, Jing & Borner 1997; Ma & Fry 2000; Seljak 2000), the halo occupation 
models (e.g., Jing, Mo & Borner 1998; Peacock & Smith 2000; Zheng et al. 2005), and the 
conditional luminosity function (e.g., Yang, Mo & van den Bosch 2003). 

Recently Gao et al (2005), using a very large iV-body simulation of a ACDM cosmology, 
reexamined the dependence of halo bias on the properties of dark matter halos. They found 
that, for low- mass halos at redshift z = with <C M* 3, the bias depends not only 
on the mass but also on the formation time of dark matter halos. For a given mass, halos 
that have earlier formation times are on average more strongly clustered in space. This age- 
dependence of halo clustering has been confirmed by a number of independent studies (e.g. 
Zhu et al. 2006; Harker et al. 2006). Since halo formation time is known to be correlated 
with halo concentration (e.g. Navarro, Frenk & White 1997; Jing 2000; Wechsler et al 2002; 
Zhao et al. 2003a,b; Lu et al. 2006), the age dependence of halo clustering also shows up as 
a halo-concentration dependence of halo clustering (Wechsler et al. 2005). For massive halos 
with M n 3> M*, the age- dependence may be reversed, as shown in Wechsler et al. (2005) 
and Wetzel et al. (2006). 

Although the mass-dependence of halo clustering is well understood in the excursion- 
set model based on spherical collapse in Gaussian initial density field (Mo & White 1996), 
the formation-time dependence (or age-dependence) is not. In the simplest excursion-set 
model from which the Press-Schechter formula is derived (Bond et al. 1991), the formation 
history of a halo is entirely determined by the local density field, independent of the mass 
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1 M t is the characteristic mass scale at which the RMS of the linear density field is equal to 1.686 at the present time. For the 
present simulation M* 1.0 X 10 13 h~ 1 MQ. 
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distribution on larger scales. Therefore, the large-scale clustering of halos is independent 
of formation time in this model. An age-dependence of halo clustering is expected in the 
excursion-set model based on ellipsoidal collapse (Sheth, Mo & Tormen 2001), because here 
the collapse of a halo depends not only on the local density field but also on the tidal field 
generated by its large-scale environment. Unfortunately, this aspect of the model was not 
explored in any detail in Sheth et al. (2001). Instead, these authors used a mean relation 
between halo mass and the shape parameter of the large-scale tidal field to derive an average 
critical collapse overdensity that depends only on halo mass. 

Since the properties of galaxies that form in a halo are expected to depend on the forma- 
tion history of the halo, understanding the origin of the age-dependence of halo clustering is 
crucial for a better understanding about how galaxies are related to their large-scale environ- 
ments. In this paper we use a high-resolution iV-body simulation to study how the formation 
of dark matter halos is affected by their environments, and how such environmental effects 
can be used to understand the age-dependence of halo clustering. The outline of the paper 
is as follows. In Section [2] we describe briefly the simulation to be used and how dark halos 
are identified. In Section [3] we study the age-dependence of halo clustering in terms of halo 
properties in the initial density field. In Section H] we examine the environments of dark 
halos and how they may affect the formation histories of dark matter halos. We discuss and 
summarize our results in Section [51 

2 SIMULATION AND DARK MATTER HALOS 

The simulation used in this paper was obtained using the P 3 M code described in Jing & 
Suto (2002). This simulation assumes a spatially-flat ACDM model, with density parameters 
Q m = 0.3 and = 0.7, and the cold dark matter power spectrum given by Bardeen et al 
(1986), with a shape parameter T = Q m h = 0.2 and an amplitude specified by as = 0.9. 
The CDM density field was represented by 512 3 particles, each having a mass of M p ~ 
6.2 x 10 8 M Q , in a cubic box of 100 /i~ 1 Mpc. The softening length is 10/i _1 kpc. Outputs 
are made at 60 different redshifts between z = 15 and z = 0, with an interval given by 
Aln(l + z) = 0.047. 

Halos are identified with the friends-of-friends algorithm (e.g. Davis et al. 1985) with a 
linking length 0.2 times the mean particle separation. In order to ensure that only physically 
bound systems are selected, we discard all particles selected into a halo but are not bound on 
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the basis of their binding energies. However, this step does not have any significant impact 
on our results. Halos at z = are related to their progenitors at higher z through halo 
merging trees. A halo in an earlier output is considered to be a progenitor of the halo in 
question if more than half of its particles are found in the final halo. The formation time of 
a halo is defined as the time when its most massive progenitor has exactly half of the final 
mass. Interpolations between adjacent outputs were adopted when estimating the formation 
times. In Fig. [1] we show the distribution of formation redshift for halos in two mass ranges. 
These results are similar to those obtained in Li et al. (2005). 



3 THE AGE DEPENDENCE OF HALO CLUSTERING 

In order to examine the age-dependence of halo clustering in the present simulation, we 
estimate the mean overdensity of dark matter within a sphere of a given radius R around 
each dark matter halo, Sh{R) = M(R)/M(R) — 1, where M(R) is the mass within a radius 
R around the halo, and M(R) = (4nR 3 /3)]) with p being the mean mass density of the 
universe. We estimate the bias parameter of a given set of halos using 

(««)> ( ' 

where (6h(R)) is the average overdensity over the set of halos in consideration, and (6 m (R)) is 
the average overdensity within all spheres of radius R centered on dark matter particles. This 
definition of b is different from that in Gao et al. (2005) where the auto-correlation function 
of halos is used to derive the bias factor. Our test showed that these two definitions lead 
to similar results. The advantage of our definition is that the statistic is more stable owing 
to the fact that there are many more dark matter particles than halos in the simulation. 
Fig. [2] shows the bias parameter obtained in this way as a function of halo formation time 
for halos of different masses. Results shown are the average for five different choices of 
R: 5, 7.5, 10, 12.5 and 15 /i _1 Mpc. We see clearly that the bias factor increases with halo 
formation redshift for low-mass halos (M^ < M+) and the dependence becomes weaker for 
Mh ~ M*. These results are in good agreement with those obtained by Zhu et al. (2006), 
using the same set of simulations, from the auto-correlation functions of dark halos. For 
massive halos, there is indication that the age-dependence may be reversed. Unfortunately, 
the present simulation is too small to allow us to study the age-dependence in detail for 
massive halos with > 10 13 /i~ 1 M Q . Note that part of the inverse relation between halo 
bias and formation time for massive halos is due to the use of non-zero mass bin. More 
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Figure 1. The distribution of formation redshift for halos in two mass ranges, 1.2 X 10 11 /i _1 M Q > M h > 6.2 X lO lo ft _1 M0 
(solid curve) and M h > 1.0 X 10 13 /i _1 M Q (dashed curve). 

massive halos within such a bin have, on average, larger bias due the the strong mass- 
dependence of halo bias for massive halos. Since the more massive halos are on average also 
younger, an inverse relation between halo bias and formation time can be produced, even if 
age- dependence of halo bias is completely absent for massive halos of a given mass. 

In order to understand the origin of the age-dependence of halo clustering, we examine 
the properties of dark matter halos in the initial density field. For each halo identified at the 
present time [z — 0), we trace the positions of all its particles back to the initial condition, 
and define the initial center of the halo to be the barycenter of these particles in the initial 
density field (set at z\ = 72). At each center, we find the largest spherical region within 
which the initial average overdensity exceeds the critical value 5 SC for spherical collapse at 
z = 0. We call this sphere the initial patch of the halo, and denote the mass within it by M ;1 . 
For a given halo, its initial center may be contained in the initial patches of other halos, and 
we denote the mass within the largest patch by M±. By definition M ;1 < M±. According to 
Press-Schechter theory (Press & Schechter 1974), the mass of the halo in which the center 
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Figure 2. The bias factor, defined by Eq. [T] as a function of halo formation redshift. Results are shown for halos in different 
mass ranges, as indicated in the panel. 



in consideration will eventually end up should be equal to Mi instead of M ;1 . The reason 
for this is discussed in detail in Sheth, Mo & Tormen (2001). In Figs. [3] we show the ratio 
between the initial mass Mi and the final halo mass, Mh, versus halo formation time for 
halos in two mass ranges. These plots show that, for many low- mass halos, the initial mass 
assigned according to the spherical collapse model is much larger than the final halo mass. 
Note that the initial/final mass ratio can be as large as 1,000 for some low-mass halos. 
This discrepancy between the Press-Schechter theory and cosmological simulations has been 
discussed in considerable detail in Sheth et al. (2001), and is the primary reason why Press- 



Environmental Dependence of Cold Dark Matter Halo Formation 7 





2 4 6 

Formation Redshift 



0.0 



0.5 1.0 1.5 2.0 
Formation Redshift 



2.5 



Figure 3. The ratio between initial mass, Mj, and the final mass as a function of halo formation redshift. Left panel is for 
low-mass halos with 1.2 X 10 11 /i _1 M Q > M h > 6.2 X lO lo /l -1 M , while right panel is for halos with M h > 1.0 X 1O 13 /i _1 M . 
The three curves in each panel represent the median, 20 and 80 percentiles. 

Schechter formula over-predicts the abundance of low-mass halos. For high-mass halos, the 
ratios are much closer to 1, suggesting that Press-Schechter formalism works better for 
high-mass halos. 

Fig.[3]also shows a clear trend that the initial/final mass ratio on average increases with 
the formation redshift, and the trend is stronger for halos of lower masses. According to the 
halo bias model of Mo & White (1996), which is based on identifying halos in the initial 
density field according to spherical collapse model, halos of larger masses should be more 
strongly correlated. It is therefore interesting to check whether the age-dependence of halo 
clustering observed in the simulation can be explained by the correlation between initial 
mass and formation time obtained here. To do this, we select low-mass halos with final 
masses 1.2 x 10 11 /i~ 1 M Q > Mh > 6.2 x 10 10 /i -1 M Q , and massive halos with final masses 
Mh > 10 13 /i _1 M Q . We bin halos in each of the two mass ranges according to their initial 
masses Mi, and compute the bias factor as a function of M;. The results are shown in 
Fig.Hl As expected, the bias factor increases strongly with mass. The prediction of the Mo 
& White (1996) bias model is shown as the solid curve. The model matches the simulation 
result for M; 3> M*, suggesting that large spherical patches are strongly clustered regardless 
of whether or not a small halo will eventually form in them. However, for M- x ^ M+, the 
model overpredicts the bias factor. This is not surprising, because small patches in which a 
low-mass halo form at the present time are biased towards low-density regions. 

We can assign to each halo a bias factor, b, according to its initial mass M 4 using the bias- 
mass relation given by the simulation. We can then average the values of b for halos in bins 
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Figure 4. Bias factor as a function of initial mass Mi for low-mass halos with 1.2 X 10 n ft _1 M s > M h > 6.2 X 10 10 /i _1 M Q 
(squares), and massive halos with M h > 1.0 x IO^H^Uq (t riangles). The solid curve shows the prediction of the Mo & White 
(1996) bias model. 

of formation redshift to predict the age-dependence of halo clustering. The dashed curves in 
Fig. [5] show the bias-formation redshift relations obtained in this way for halos in two mass 
ranges. The predictions match well the simulation results, suggesting that the observed age- 
dependence of clustering for halos of a given final mass can be entirely explained by the 
difference in the initial mass. Old, low-mass halos are more strongly clustered than their 
younger counterparts simply because they are associated with perturbations of higher mass 
(i.e. patches corresponding to higher peaks) in the initial density field. The age-dependence is 
weaker for more massive halos, because their M; represent their final masses more faithfully. 
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Figure 5. Bias factor as a function of formation redshift. Symbols connected with solid lines are the results obtained from 
simulations (the same as those shown in Fig. [2), while symbols connected by dashed lines are the corresponding predictions 
based on the correlation between bias and initial mass (see text for detail). Results are shown for low-mass halos with 1.2 X 
1O 11 /i~ 1 M > M h > 6.2 X lO lo /i _1 M , and for massive halos with M h > 1.0 X 1O 13 /i _1 M , as indicated in the panel. For 
the other two mass ranges shown in Fig 2, the model predictions match the simulation results as well. They are not plotted 
here for the sake of clarity. 



4 THE ENVIRONMENTAL EFFECTS OF HALO FORMATION 

The question is, of course, why some low-mass halos can survive in a patch which itself should 
collapse to form a virialized halo according to spherical collapse model, and what prevents 
such halos from accreting mass in later times. To have some ideas about what is going on 
during the formation of a low-mass old halo, we show in Fig.[6] the spatial distribution of 
dark matter particles in the neighborhood of such a halo for a number of representative 
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snapshots, with the particles that eventually end up in the final halo marked in colour. As 
one can see, the main part of the small halo collapses quickly to form a halo, which then 
grows only slowly while orbiting around a bigger halo that forms later. The early collapse 
is driven largely by the relatively high initial density of the material associated with the 
halo, as is demonstrated in Fig. [7] which shows the initial overdensity around a halo center 
as a function of halo formation redshift. In Fig.[7J the overdensity is defined as the mean 
within a sphere that is centered at the initial position of the halo and has a radius such 
that the total mass it contains is equal to the mass of the halo. For low-mass halos, there is 
a strong, positive correlation between the mean overdensity and the formation time. Halos 
with earlier formation times have initial overdensity much higher than the critical density 
for spherical collapse (indicated as the horizontal line in the figure). With such overdensities, 
these perturbations can collapse in early times to form virialized halos. The situation is very 
different for massive halos, where the mean overdensity are all close to the critical value for 
spherical collapse [see Fig. [7(b)]. 

There are two possibilities why low-mass old halos stop accreting mass at relatively early 
times. The first is that these halos reside in locally low-density regions, and so there is no 
material for them to accrete. To test this, we estimate the total mass between r vir and 3r V j r 
(with r v ; r the virial radius) for each halo and plot the mean mass density (in units of mean 
density of the universe) within these two radii versus halo formation redshift in Fig.[8](a). 
This figure shows that there is in fact a slightly larger amount of mass around older halos, 
and so the truncation of the growth of old, low-mass halos cannot be explained by the lack 
of material around them. However, the availability of material is a necessary, not a sufficient 
condition for halo accretion. In order for the material around a halo to be accreted by the 
halo, the velocity of the material to be accreted relative to the halo must be sufficiently low. 
To see if this condition is fulfilled, we estimate the velocity difference |v p — Vh| for all particles 
with halo-centric distances between r v - ir and 3r v i r , where Vh is the velocity of the host halo, 
and v p is the velocity of a CDM particle. Fig.[8](b) shows the velocity difference (in units of 
halo virial velocity) versus halo formation redshift. There is a clear trend that particles in 
the neighborhood of an older halo have systematically larger velocity differences with the 
host halo. For the oldest halos, the velocity differences for most particles are significantly 
larger than the virial velocities of the halos. These particles cannot be accreted by the halo, 
because they are too energetic. Thus, the reason why late accretions by old, low-mass halos 
are truncated is that they are embedded in 'hot' regions where particles are moving too 



Environmental Dependence of Cold Dark Matter Halo Formation 11 

fast for them to capture. The situation is very different for massive halos [Figs.lHfc) and 
(d)]. First of all, the mean densities in their neighborhoods are quite independent of halo 
formation redshift [Fig.[S(c)]. More importantly, the velocity differences between a halo and 
particles around it are smaller than the virial velocity and so most of these particles can be 
accreted by the halo in subsequent evolution. 

It is easy to understand why old, low-mass halos are embedded in hot environments. As 
shown in Fig. [HI these halos are formed in the vicinity of big structures that are dominated 
by one or more massive halos. To see this more clearly, we plot in Fig. [9] the mass of the 
most massive neighboring halo within a distance of 2/i _1 Mpc from a low-mass halo as a 
function of the formation redshift of the low-mass halo. An old low-mass halo almost always 
has a massive halo in its neighborhood. These old low-mass halos are similar to the subhalo 
population see in high-resolution A-body simulations (e.g. Klypin et al. 1999; Moore et 
al. 1999), except that they have not yet merged into the virialized part of a massive halo. 
Although the formation of the more massive neighboring halos occurred at later times, the 
tidal field of the large structure can accelerate the particles around low-mass halos, thereby 
increasing their velocities relative to the low-mass halos. Furthermore, massive pancakes and 
filaments can form in high-density patches prior to the formation of massive halos, and heat 
up the particles in them in a way similar to the preheating mechanism discussed in Mo et 
al. (2005). Note that these two processes are related, because the formation of pancakes and 
filaments is related to the local tidal field. As discussed in Mo et al. (2005), the tidal effect 
and preheating are expected to have more significant impact on smaller halos, which may 
explain why the age-dependence of halo clustering is weaker for more massive halos. 

5 DISCUSSION AND SUMMARY 

In this paper, we have used a high-resolution A-body simulation to study how the formation 
of CDM halos may be affected by their environments, and how such environmental effects 
may produce the age- dependence of halo clustering seen in recent large A-body simulations. 
We have shown that old low-mass halos are almost always associated with initial perturba- 
tions that correspond to higher peaks than their present masses imply. Consequently, these 
halos are strongly clustered in space despite of their low masses at the present time. Because 
of their associations with large density perturbations, old low-mass halos are almost always 
found in the vicinity of big structures that are dominated by one or more massive halos at 
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the present time. Mass accretions into these low-mass halos are limited at late times by the 
tidal field of the larger neighboring structures. Such environmental effects are weaker for 
more massive halos, which explains why the age-dependence of halo clustering is weaker for 
more massive halos. 

Our results demonstrate clearly that environmental effects play an important role during 
the formation of dark matter halos, especially of low-mass halos. Since the properties of 
galaxies that form in a halo are expected to depend on the formation history of the halo, 
such environmental effects may have important implications for galaxy formation in different 
environments. Galaxies that form in old, low- mass halos are presumably old and faint, and 
our results suggest that many such galaxies may be located in the vicinity of relatively 
massive galaxy systems. These galaxies likely represent an extension of the satellite galaxies 
observed in galaxy systems such as the local group and other groups/clusters of galaxies. 
Observationally, it is known that satellite galaxies in galaxy groups are dynamically old 
systems (dwarf ellipticals and dwarf spheroids) that contain mainly old stellar populations. 
It is interesting to see whether the stellar populations and kinematics of low-mass galaxies in 
the neighborhoods of nearby galaxy groups have properties similar to those of the satellite 
population. 

Theoretically, it is very interesting to see if the ellipsoidal collapse model that incorpo- 
rates large-scale tidal field into the dynamics can indeed explain the environmental effects 
discussed here and the age-dependence of halo clustering observed in iV-body simulations. 
We will come back to this problem in a future work. 
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Figure 6. CDM particle distributions around a low-mass old halo (with final mass M% = 1.2 X 10 1 h~ Mg and formation 
redshift 3.86) at a number of different redshifts (as indicated in the panels). Particles that eventually end up in the final halo 
are plotted in red. 




Figure 7. The correlation between the initial overdensity, defined on the halo mass scale, and halo formation redshift. Results 
are shown for low-mass halos with 1.2 X lO 11 fe- 1 M > M h > 6.2 X IO^H^Mq (left panel) and for massive halos with 
M h > 1.0 X 10 13 /i _1 M Q (right panel). The three curves in each panel show the median, 20 and 80 percentiles. The horizontal 
lines indicates the critical overdensity for collapse (at the initial redshift z = 72) based on spherical collapse model. 
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Figure 8. Left panels: the mean density (in units of the mean density of the universe) between r v i r and 3r v i r as a function of 
halo formation redshift. Each point represents a halo. Right panels: the velocity difference for particles located between r v i r and 
3r v i r as a unction of halo formation redshift. Here each point represents a dark matter particle. Upper panels are for low-mass 
halos with 1.2 X lO n /i _1 M > M h > 6.2 X 1O 1o /i _1 M , and lower panels and for massive halos with M h > 1.0 X 10 13 fe _1 M Q . 
The three curves in each panel are the median, 20 and 80 percentiles. 
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Figure 9. The ratio between Mi argest and as a function of halo formation redshift. Here Mi argcst is defined to be the mass 
of the most massive halo that is within 2/i -1 Mpc from the halo in consideration. The three curves show the median, 20 and 
80 percentiles. 



